Bianchi type I universe with viscous fluid 



o 
O 

m 



Bijan Saha 

Laboratory of Information Technologies 
Joint Institute for Nuclear Research, Duhna 
14 1980 Duhna, Moscow region, Russi^ 
(Dated: February 4, 2008) 

We study the evolution of a homogeneous, anisotropic Universe given by a Bianchi 
type-I cosmological model filled with viscous fluid, in the presence of a cosmological 
' constant A. The role of viscous fluid and A term in the evolution the BI space-time 

^ ■ is studied. Though the viscosity cannot remove the cosmological singularity, it plays 

, a crucial part in the formation of a qualitatively new behavior of the solutions near 

singularity. It is shown that the introduction of the A term can be handy in the 
elimination of the cosmological singularity. In particular, in case of a bulk viscosity, 
it provides an everlasting process of evolution (A < 0), whereas, for some positive 
values of A and the bulk viscosity being inverse proportional to the expansion, the 
BI Universe admits a singularity-free oscillatory mode of expansion. In case of a 
' constant bulk viscosity and share viscosity being proportional to expansion, the 

■ model allows oscillatory mode accompanied by an exponential growth even with a 
. negative A. Space-time singularity in this case occurs at t — > —00. 

^ . PACS numbers: 03.65.Pm, 04.20.Jb, and 04.20.Ha 

■r^lj- ' Keywords: Bianchi type I (BI) model, Cosmological constant, viscous fluid 

o ■ 

I. INTRODUCTION 

^ ' The investigation of relativistic cosmological models usually has the energy momentum 

■ tensor of matter generated by a perfect fluid. To consider more realistic models one must take 
into account the viscosity mechanisms, which have already attracted the attention of many 
researchers. Misner [l[ |^ suggested that strong dissipative due to the neutrino viscosity 
may considerably reduce the anisotropy of the blackbody radiation. Viscosity mechanism 
in cosmology can explain the anomalously high entropy per baryon in the present universe 
[1, i^J . Bulk viscosity associated with the grand- unified-theory phase transition may lead 
to an inflationary scenario 0, 0, S] • 

A uniform cosmological model filled with fluid which possesses pressure and second (bulk) 
viscosity was developed by Murphy [9]. The solutions that he found exhibit an interesting 
feature that the big bang type singularity appears in the infinite past. Exact solutions of 
the isotropic homogeneous cosmology for open, closed and flat universe have been found by 
Santos et al 10], with the bulk viscosity being a power function of energy density. 

The nature of cosmological solutions for homogeneous Bianchi type I (BI) model was 



investigated by Belinsky and Khalatnikov [11| by taking into account dissipative process 



due to viscosity. They showed that viscosity cannot remove the cosmological singularity 
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but results in a qualitatively new behavior of the solutions near singularity. They found 
the remarkable property that during the time of the big bang matter is created by the 
gravitational field. BI solutions in case of stiff matter with a shear viscosity being the power 
function of of energy density were obtained by Banerjee whereas BI models with bulk 
viscosity (rj) that is a power function of energy density e and when the universe is filled 
with stiff matter were studied by Huang The effect of bulk viscosity, with a time 

varying bulk viscous coefficient, on the evolution of isotropic FRW models was investigated 
in the context of open thermodynamics system was studied by Desikan jl^]. This study was 
further developed by Krori and Mukherjee for anisotropic Bianchi models. Cosmological 
solutions with nonlinear bulk viscosity were obtained in Models with both shear and 
bulk viscosity were investigated in |l7L Il8l|. 

Though Murphy Q claimed that the introduction of bulk viscosity can avoid the initial 
singularity at finite past, results obtained in show that, it is, in general, not valid, since 
for some cases big bang singularity occurs in finite past. 

We studied a self-consistent system of the nonlinear spinor and/or scalar fields in a BI 
spacetime in presence of a perfect fluid and a A term [20, |2l| in order to clarify whether the 
presence of a singular point is an inherent property of the relativistic cosmological models 
or it is only a consequence of specific simplifying assumptions underlying these models. 
Recently we have considered a nonlinear spinor field in a BI Universe filled with viscous 
fluid [22]. Since the viscous fluid itself presents a growing interest, we study the influence of 
viscous fluid and A term in the evolution of the BI Universe in this report. 



II. DERIVATION OF BASIC EQUATIONS 

Using the variational principle in this section we derive the fundamental equations for 
the gravitational field from the action ()2.1|) : 

S{g;e) = J C^gdQ (2.1) 

with 

C -Cgrav. ~l~ -^vf- (2-2) 

The gravitational part of the Lagrangian (j2.2j) £grav. is given by a Bianchi type-I metric, 
whereas the term £vf describes a viscous fluid. 

We also write the expressions for the metric functions explicitly in terms of the vol- 
ume scale r defined bellow fl2.2Hj) . Defining Hubble constant ()2.H7j) in analogy with a flat 
Friedmann-Robertson- Walker (FRW) Universe, we also derive the system of equations for 
r, H and e, with e being the energy density of the viscous fluid, which plays the central role 
here. 



A. The gravitational field 

As a gravitational field we consider the Bianchi type I (BI) cosmological model. It is 
the simplest model of anisotropic universe that describes a homogeneous and spatially flat 
space-time and if filled with perfect fluid with the equation of state p = (e, ( < 1, it 
eventually evolves into a FRW universe [2^. The isotropy of present-day universe makes 
BI model a prime candidate for studying the possible effects of an anisotropy in the early 
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universe on modern-day data observations. In view of what has been mentioned above we 
choose the gravitational part of the Lagrangian ()2.2j) in the form 



(2.3) 



where R is the scalar curvature, k = SnG being the Einstein's gravitational constant. The 
gravitational field in our case is given by a Bianchi type I (BI) metric 



ds"^ = dt^ — a^dx^ — iP'dy^ — c^dz^ 



(2.4) 



with a, 6, c being the functions of time t only. Here the speed of light is taken to be unity. 
The metric ()2.4|1 has the following non-trivial Christoffel symbols 
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The nontrivial components of the Ricci tensors are 
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From ()2.(i|) one finds the following Ricci scalar for the BI universe 

R 

The non-trivial components of Riemann tensors in this case read 
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(2.5) 



(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 



(2.7) 



(2.8) 



Now having all the non-trivial components of Ricci and Riemann tensors, one can easily 
write the invariants of gravitational field which we need to study the space-time singularity. 
We return to this study at the end of this section. 
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B. Viscous fluid 



The influence of tlie viscous fluid in the evolution of the Universe is performed by means 
of its energy momentum tensor, which acts as the source of the corresponding gravitational 
field. The reason for writing C^f in ()2.2|) is to underline that we are dealing with a self- 
consistent system. The energy momentum tensor of a viscous field has the form 

^/Tm = + P')u^u'' - + -qg^^lu^-fs + Up.,^ - u^u^Up-a - Upu'^u^.a], (2.9) 

where 

p' = P-(e-^^K- (2.10) 

Here e is the energy density, p - pressure, t] and ^ are the coefficients of shear and bulk 
viscosity, respectively. Note that the bulk and shear viscosities, rj and ^, are both positively 
definite, i.e., 

r/>0, ^>0. (2.11) 
They may be either constant or function of time or energy, such as : 

r] = \A\e'', ^=\B\e^. (2.12) 

The pressure p is connected to the energy density by means of a equation of state. In this 
report we consider the one describing a perfect fluid : 

p = Ce, Ce(0,l]. (2.13) 

Note that here ( 0, since for dust pressure, hence temperature is zero, that results in 
vanishing viscosity. 

In a comoving system of reference such that m'^ = (1, 0, 0, 0) we have 

^oV) = ^' (2-14a) 
T^M = -P' + ^V^, (2.14b) 

T^Im = -P' + ^V^ (2.14c) 

T^Im = -p' + 2v-^. (2.14d) 

Let us introduce the dynamical scalars such as the expansion and the shear scalar as 
usual 

9 = u^^, a2 = la^,a^-, (2.15) 

where 



+ - loP,,. (2.16) 

Here P is the projection operator obeying 

P^ = P. (2.17) 
For the space-time with signature (+, — , — , — ) it has the form 

P^u = 9f..-u^u,, Pi: = 5^, - u^u,. (2.18) 
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For the BI metric the dynamical scalar has the form 



and 



+ ^ + H = I (2.19) 
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C. Field equations and their solutions 

Variation of ()2.H) with respect to metric tensor g^^ gives the Einstein's field equation. In 
account of the A-term we then have 

G'^i, = R^, ^^1^ 
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rewrite the Eq. (|2.21|) as 




= kTI~ 


A, 


(2.22a) 


= KTi - 


A, 


(2.22b) 


= - 


A, 


(2.22c) 


= - 


A, 


(2.22d) 



where over dot means differentiation with respect to t and T^^ is the energy-momentum 
tensor of a viscous fluid given above ()2.14|) . 



1. Expressions for the metric functions 

To write the metric functions explicitly, we define a new time dependent function r(t) 

r = abc = (2.23) 

which is indeed the volume scale of the BI space-time. 

Let us now solve the Einstein equations. In account of ()2.14p subtracting ()2.22ap from 
fl2.22b|) . one finds the following relation between a and b 

- = DiexpiXi dtj. (2.24) 

Analogically, we find 

-=D2expiX2 dtl - = D3expiXs dt\. (2.25) 
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Here Di, D2, D^, Xi, X2, X3 are integration constants, obeying 

D1D2D3 = 1, Xi + X2 + X3 = 0. (2.26) 
In view of ()2.26|1 from ()2.24|1 and ()2.25|) we write the metric functions explicitly [l^l 



a{t) 
b{t) 
c{t) 



AiT^^^ exp 
A2T^^^ exp 
A^T^^^ exp 
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where 



A^ = ^(D^/Ds), A2 = yi/iDlDs) 
Bi = Xi— X3, B2 = —(2X1 + X3), 



Pi/^i), 
-B3 = Xi + 2X3. 



(2.27a) 
(2.27b) 
(2.27c) 



Thus, the metric functions are found explicitly in terms of r and viscosity. 

As one sees from ()2.27ap . ()2.27b|) and ()2.27cp . for r = with n > 1 the exponent tends 
to unity at large t, and the anisotropic model becomes isotropic one. 



2. Singularity analysis 



Let us now investigate the existence of singularity (singular point) of the gravitational 
case, which can be done by investigating the invariant characteristics of the space-time. In 
general relativity these invariants are composed from the curvature tensor and the metric 
one. In a 4D Riemann space-time there are 14 independent invariants. Instead of analyzing 
all 14 invariants, one can confine this study only in 3, namely the scalar curvature Ji = i?, 
I2 = R/.uR^", and the Kretschmann scalar h = RaPt^uR'"'^^'' Qll^. At any regular space- 
time point, these three invariants /i, I2, I3 should be finite. Let us rewrite these invariants 
in detail. 

For the Bianchi I metric one finds the scalar curvature 



h = R 



d b c db be cd 
- + - + - + -- + -- + — 
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(2.28) 



Since the Ricci tensor for the BI metric is diagonal, the invariant I2 = R^^R^^ 
sum of squares of diagonal components of Ricci tensor, i.e.. 



RI^Ri^ is a 



with the components of the Ricci tensor being given by ()2.6p . 

Analogically, for the Kretschmann scalar in this case we have J3 = R^'^ ^ 
squared components of all nontrivial R^^ which in view of ()2.8|1 can be written as 



(2.29) 



R'\,y+{R'' 
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Let us now express the foregoing invariants in terms of r. From Eqs. ()2.27|) we have 

r g^2K J r]dt 

ai = Ay^^expi (5,/3) / ^^dt ] , (2.31a) 



r{t) 

- = k (2 = 1,2,3,), (2.31b) 

Oj or 

di _ Srf -2t^ -rBte-^^^'T^^ -6Kr]rBie-'^''^'T^^ + Bfe-^^^'^'^^ 



(2.31c) 



i.e., the metric functions a,b,c and their derivatives are in functional dependence with r. 
From Eqs. ()2.3ip one can easily verify that 



h oc \, hoc ^, I3 oc ^. 

^4 ^4 



Thus we see that at any space-time point, where r = the invariants Ji, I2, and Is become 
infinity, hence the space-time becomes singular at this point. 



D. Equations for determining r 

In the foregoing subsection we wrote the corresponding metric functions in terms of 
volume scale r. In what follows, we write the equation for r and study it in details. 
Summation of Einstein equations (I2.22a|) . (I2.22bjl . (I2.22c|l and 3 times ()2.22djl gives 

3 3 

f--K^T = -K{e-p)T-3AT. (2.32) 

For the right-hand-side of ()2.32|) to be a function of r only, the solution to this equation is 
well-known [26j. 

Let us demand the energy- momentum to be conserved, i.e., 

TH;. = t;^u + r;:,T; - f^^t; = o, (2.33) 

which in our case has the form 

i(rT°)-^T^-^T|-^r| = 0. (2.34) 
After a little manipulation from ()2.34|) we obtain 



£ + - (e + -r/)^ + H>^T^ - A) = 0, (2.35) 



T 3 

where 

LJ = e + p, (2.36) 

is the thermal function. 

Let us now, in analogy with Hubble constant in a FRW Universe, introduce a generalized 
Hubble constant H : 

f a h c , , 

- = - + - + - = 3H. 2.37 

T a b c 
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Then ()2.32|) and ()2.35p in account of ()2.14|) can be rewritten as 

ij = I [3^H -uj)- {3H^ -KE + A), (2.38a) 

e = 3H{3^H -io) +4r]{3H^ - Ke + A). (2.38b) 

The Eqs. fl2.38|) can be written in terms of dynamical scalar as well. 
In account of ()2.27p one can also rewrite share scalar ()2.20|) as 

2a' = 6(-^i+-^i^3 + X|) ^_4.j,.._ (2.39) 

From ()2.22djl now yields 

1^2 _ ^2 ^ _ ^ (2.40) 
The Eqs. ()2.38|) now can be written in terms of 6 and a as follows 



e = ^[^e~u) -3a\ (2.41a) 
e = e{ie-u) +4r/(Tl (2.41b) 

Note that the Eqs. ()2.41|) coincide with the ones given in 3 . 



E. Some special solutions 

In this subsection we simultaneously solve the system of equations for r, H, and e. It is 
convenient to rewrite the Eqs. ()2.37j) and ()2.38j) as a single system : 

f = 3Ht, (2.42a) 

if = I [3iH -uo)- (3iJ^ -KE + A), (2.42b) 

e = 3H{3^H -uj)+ At]{3H' -ke + A). (2.42c) 

In account of ()2.36|) . ()2.12|) and ()2.13p the Eqs. ()2.42|) now can be rewritten as 

f = 3Ht, (2.43a) 

if = ^ [3Be'^H - (1 + C)^) - (3//2 -Ks + A), (2.43b) 

e = 3H[3Be^H - (1 + C)^) + 4Ae° (3/7^ -ke + A). (2.43c) 

The svster n (I2.42|) have been extensively studied in hterature either partially [T^, 
or in general [ll| . In what follows, we consider the system ()2.42j) for some special choices of 
the parameters. 

1. Case with hulk viscosity 



Let us first consider the case when the real fluid possesses the bulk viscosity only. The 
corresponding system of Eqs. can then be obtained by setting r/ = in (j2.42j) or A = in 



Bianchi type I universe with viscous Ruid 



9 



()2.43|) . In this case the Eqs. (I2.42a|l and (I2.42bjl remain unaltered, while fl2.42c|) takes the 
form 

e = 3H{3^H -uj). (2.44) 
In view of ()2.44|) the system ()2.42j) admits the following first integral 

{ke - 3H^ - A) = Ci, Ci = const. (2.45) 

The relation ()2.45|) can be interpreted as follows. At the initial stage of evolution the volume 
scale T tends to zero, while, the energy density e tends to infinity. Since the Hubble constant 
and the A term are finite, the relation ()2.45|) is in correspondence with the current line of 
thinking. Let us see what happens as the Universe expands. It is well known that with 
the expansion of the Universe, i.e., with the increase of r, the energy density e decreases. 
Suppose at some stage of expansion r — oo and e — > 0. Then from ()2.45|) follows that at 
the stage in question 

3H^ + A ^ 0. (2.46) 

In case of A = 0, we find H = 0, i.e., in absence of a A term, once r —>■ oo, the process of 
evolution is terminated. As one sees from ()2.4fi|l . for the H to make any sense, the A term 
should be negative. In presence of a negative A term the evolution process of the Universe 
never comes to a h alt, it either expands further or begin to contract depending on the sign 
of H = ±^-A/3, A < 0. Thus we see that the Universe may be infinitely large only if 
A < 0. 

Let us now consider the case when the bulk viscosity is inverse proportional to expansion, 
i.e., 

= C2, C2 = const. (2.47) 

Now keeping into mind that 9 = f/r = 3H, also the relations ()2.42aj) . ()2.3(jj) and ()2.13|1 the 
Eq. ()2.44j) can be written as 

From the Eq. ()2.48j) one finds 



^ (2.48) 



C2-(l + C)^ r 

[C2 + CsT-'-'^^^] , (2.49) 



with C3 being some arbitrary constant. Further, inserting e from ()2.49|1 into ()2.32|1 one finds 
the expression for r explicitly. 

Taking into account the equation of state ()2.13|) in view of ()2.47|) and ()2.49|) . the Eq. 
()2.32|1 admits the following solution in quadrature : 



/ 



dr 

t + to, (2.50) 



where C| and to are some constants. Further we set to = 0. Here, Cq = 3kC2/{1 + C) — 3A 
and Cl = 3kC3/{1 + (). As one sees, Cq is negative for 

A>kC2/{1 + C)- (2.51) 

It means that for a positive A obeying ()2.51|) (we assume that the constant C2 is a positive 
quantity) r should be bound from above as well. Let us now rewrite the Eq. ()2.50|1 in the 
form 

f = V2[^-W(r)], (2.52) 
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where E = C|/2 can be viewed as energy and i/(r) = —0.5{CqT^ +CIt^~'^) can be seen as the 
potential [cf. Fig. P corresponding to the Eqn. ()2.32j) . Depending on the value of E there 
exists two types of solutions: for E > we have non-periodic solutions, i.e., after reaching 
some maximum value (say Tmax) the BI Universe begin to contract and finally coUops into 
a point, thus giving rise to a space-time singularity; for E < Bl space-time admits a 
singularity- free oscillatory mode of expansion [cf. Fig. EJ. A comprehensive description 
concerning potential can be found in [21]. 




FIG. 1: View of the potential U{t). Here we FIG. 2: Evolution of the BI space-time cor- 
set «;= 1, C2 = 1, and C3 = 1. Perfect fluid responding to the potential given in Fig. 
corresponds to a radiation, i.e., ( = 0.33, The initial value of the volume scale in this 
and the Cosmological constant is taken to case is taken to be tq = 0.2. 
be A = 0.8 



As a second example we consider the case, when C = 1. From ()2.50|) one then finds 

T{t) = (exp(^ t) - Clexpi^^ t))/(2v^), C° > 0, (2.53a) 

r{t) = {Cll^\) sin (y^ t). C° < 0. (2.53b) 

Taking into account that Cg > for any non-positive A, from ()2.53aj) one sees that, in case 
of A < the Universe may be infinitely large (there is no upper bound), which is in line 
with the conclusion made above. On the other hand, Cq may be negative only for some 
positive value of A. It was shown in QA^ that in case of a perfect fluid a positive A 
always invokes oscillations in the model, whereas, in the present model with viscous fluid, 
it is the case only when A obeys ()2.51|) . Unlike the case with radiation where BI admits 
two types of solutions, the case with stiff matter allows only one type of solutions, namely 
the non-periodic one that corresponds to > in the previous case, since now potential 
U{t) = —0.5C^t'^ has its minimum at r = 0. 

2. Case with shear and bulk viscosity 

Let us now consider the general case with the shear viscosity t] being proportional to the 
expansion, i.e., 

r](xe = 3H. (2.54) 
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We will consider the case when 



V 



2k 



In this case from ()2.42b|) and ()2.42c|) one easily find 

3H^ = KS + Ci, C4 = const. 



(2.55) 



(2.56) 



From ()2.5(ij) it follows that at the initial state of expansion, when e is large, the Hubble 
constant is also large and with the expansion of the Universe H decreases as does e. Inserting 
the relation ()2.5f)|l into the Eqs. (j2.42bj) one finds 



dH 



VAH^ + BH + C 



(2.57) 



where, A = -1.5(1 + C), B = 1.5<, and C = 0.5C4(C - 1) - A. If the bulk viscosity is 
taken to be a constant one, i.e., ^ = const., then in view of the fact that A < 0, the Hubble 
constant H admits sinusoidal form, namely, [27l | 



H 



- 4AC sin (v^ t) + B 
2A ' 



if and only if 



VB^ - 4AC = (9/4)/t2^2 _ _ ^2)^^ _ ^ > 0. 



(2.58) 



(2.59) 



As one sees the inequality ()2.59p can be attained with a negative A as well. Further, from 
()2.42aj) one finds the expression for r: 



r(t) = Cq exp 



3 f -VW^^iAC cosi 
2 ^ 



-at) + B^T^t 



Ay/^ 

A graphical view of the evolution of r is given in Fig. IHl 



}]■ 



(2.60) 




FIG. 3: Evolution of the BI Universe with a bulk viscosity and a share 
viscosity. The volume scale r in this case expands exponentially with 
a small oscillation. 



As one sees from the Fig. El the exponential mode of evolution of the BI Universe is 
accompanied by an oscillation. The space-time singularity in this case arises aX t ^ —00. 
Note that the negative A gives rise to the exponential growth while the oscillation in the 
model occurs due to the viscosity. 
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III. CONCLUSION 

The role of viscous fluid and A term in the evolution of a homogeneous, anisotropic 
Universe given by a Bianchi type-I space-time is studied. It is shown that the A term plays 
very important role in BI cosmology. In particular, in case of a bulk viscosity, it provides an 
everlasting process of evolution with A being negative, whereas, for some positive values of A 
the BI Universe admits an oscillatory mode of expansion. It is also shown that a oscillatory 
mode of expansion of the BI space-time can take place with a negative A as well, though it 
is accompanied by an exponential mode. Oscillation in this case arises due to viscosity. In 
this report we only consider some special cases those provides exact solutions. For a better 
knowledge about the evolution, it is important to perform some qualitative analysis of the 
system ()2.42p . A detailed analysis of the system in question plus some numerical solutions 
will be presented soon elsewhere. 
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